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^ ■ New measurements of short-range and long-range two-particle correlations, azimuthal 
anisotropy, and event-by-event fluctuations from the STAR experiment for ^snn = 130 
and 200 GeV Au+Au collisions are summarized. Striking evidence is presented for large, 
non-statistical fluctuations in mean transverse momentum. Descriptions of the data in 
^ ; terms of phenomenological source function models are also presented. 

V 

■ 1. Introduction 

Determining whether or not color deconfinement and chiral symmetry restoration occur 
during relativistic heavy ion collisions at RHIC requires a multi-prong effort involving ex- 
periment, phenomenology, and theory. Experimentally, measurements of inclusive spectra 
and correlation distributions for hadrons and leptons are needed to study quantum in- 
terference, final state interactions, and non-statistical fluctuations which produce correla- 
tions evident in the observed momentum distributions. Detailed studies of non-statistical 
event-by-event fluctuations |I| by way of correlation measurements are possible at RHIC, 
especially with the large acceptance STAR detector 0. To help bridge the gap between 
experiment and theory, phenomenological source function models are being used by both 
experimentalists and theorists to describe the data. 

This paper summarizes recent correlation measurements and phenomenological source 
model descriptions from the STAR experiment for soft hadron production in Au+Au colli- 
sions at y/SNN = 130 and 200 GeV at RHIC. The data include results on Hanbury Brown 
and Twiss (HBT) pion interferometry, pion phase space densities, pion-kaon correlations, 
long-range correlations in pseudorapidity space, azimuthal anisotropy, and event-by-event 
fluctuations in mean transverse momentum. Data presented here were obtained with the 
STAR Time Projection Chamber (TPC) using the Level trigger detectors CTB (Cen- 
tral Trigger Barrel) and ZDC (Zero-Degree Calorimeter) M and should be regarded as 
preliminary. 
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2. Source Function Phenomenology 

Models of the space-time and momentum distribution of the hadronic freeze-out region 
in relativistic heavy ion collisions typically have the following form (for bosons) |fl: 

m T coshfy — 77) f / p ■ u(x) — a\ „ 1 1 , 
S(x,p) = T £ " jexp r y H - l| F(r,<^,r), (1) 

where x [with components r, <p, rj (source rapidity) and r] is the space-time coordinate, p 
(with components m?, y, 4>) is the on-shell momentum coordinate, my is the transverse 
mass, u(x) represents the collective velocity field which accounts for longitudinal (Bjorken) 
and transverse expansions, \x is the chemical potential, T is the local temperature, and 
F is an empirical modulation function. The latter may depend on the radial, azimuthal, 
rapidity and longitudinal proper time (r = \/t 2 — z 2 ) coordinates of the source. A special 
case of the general emission function in Eq. ([!]) is referred to as the "blast-wave" model. 
This model assumes longitudinal boost invariance, finite freeze-out duration with respect 
to a fixed longitudinal proper time tq, and integration over the observed rapidity y. The 
blast-wave model has been used to describe some of the STAR data taken in the first two 
RHIC runs as discussed below. 

3. HBT Pion Interferometry 

Preliminary results for three-dimensional Pratt-Bertsch || "out-side-long" HBT anal- 
yses of the recent ^Jsnn = 200 GeV Au+Au minimum bias collision data from STAR for 
charged pions (see contribution to these proceedings by M. Lopez-Noriega) are very simi- 
lar to the STAR ]5| and PHENIX |J results at 130 GeV. The measured three-dimensional 
correlations were fitted with the standard gaussian form, 

C(q , q„ q t ) = 1 + A exp[-R 2 q 2 - R 2 s q 2 s - R 2 q 2 } (2) 

where the relative momentum vector was evaluated in the longitudinal co-moving source 
(LCMS) frame and decomposed into components qi (parallel to beam direction), q Q (par- 
allel to pair transverse momentum) and q s (orthogonal to each). Standard cuts and 
corrections for track splitting, two-track merging in the STAR TPC, and Coulomb cor- 
rections were included as in Ref. 0. The correlation lengths^] or "radii," R ,R s ,Ri, at 
t/snn — 200 GeV are similar in magnitude and have the same centrality and transverse 
mass {m T ) dependences as the results at ^/snn = 130 GeV. Ri adheres to the Makhlin- 
Sinyukov [0] dependence on itit expected for Bjorken expansion with corresponding source 
lifetimes which range from approximately 8 fm/c for peripheral collisions to 10 fm/c for 
central events. The ratio R a j R s for the 200 GeV data is approximately unity in agreement 
with results reported at 130 GeV ||. This result is not understood theoretically but can 
be described with a brief emission duration (At) and sharp radial cut-off of the source 
geometry. 

2 The conjugate lengths of the relative momentum correlations characterize the homogeneity scale of the 
emitting source. It is customary to refer to these lengths as "radii." For expanding sources additional 
analysis is required to determine the total source geometry. 
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Another important aspect of the pion emission source which can be inferred from the 
HBT results is the phase space density (PSD). Bertsch || and Ferenc et al. [[| have 
shown that the pion transverse momentum (p?) distribution together with the volume of 
the two-pion HBT correlation function can be used to determine the spatial average of the 
pion PSD, (f). This method was applied to the minimum bias y / s/v r /v = 130 GeV Au+Au 
collision data from STAR for negative pions for seven centrality bins [defined by the 
percent ranges 0-5 (most central), 5-10, 10-20, 20-30, 30-40, 40-50, and 50-80 (peripheral) 
of the total reaction cross section]. The my distribution data, shown in Fig. [j], were 
fitted with a Bose-Einstein (BE) distribution, {A/[exp(rriT/T e ff) — 1]}, with effective 
temperatures approximately 200 MeV (solid curves). The spectra in Fig. p] increase in 
magnitude monotonically with centrality. The volumes of the three-dimensional Pratt- 
Bertsch HBT correlation functions for these seven centralities were obtained by direct 
summation of the binned data. The resulting estimates of (f) for three ranges of pt 
(defined by the limits 0.125, 0.225, 0.325, and 0.45 GeV/c) are shown in Fig. g. For 
comparison, similar estimates for central Pb+Pb collisions at the SPS from experiment 
NA49 H are indicated by the solid triangles in Fig. [2|. The results show that (f) increases 
with y/SNN from SPS to RHIC and monotonically increases with centrality at RHIC 
energies. Apparently the universal pion PSD suggested by Ferenc U for AGS and SPS 
energies is not valid at RHIC energies. 

An interesting prediction follows from the source function model of Tomasik and Heinz || 
which includes explicit longitudinal expansion effects, when fitted to the negative pion itlt 
spectra in Fig. [I] and the PSD values in Fig. |2] (solid curves). The model, when extrap- 
olated to low p T , predicts large values of (f) ~ 1 for central events, suggesting that 
significant multiparticle Bose-Einstein effects may be present at low px at RHIC. 
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Figure 1. Inclusive ir~ transverse 
mass spectra in units of GeV for 
seven centralities with BE distribu- 
tion model fits (solid curves). 
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Figure 2. Average 7r~ PSD versus px for seven cen- 
tralities from STAR and from NA49 for top central 
events (solid triangles) with BE distribution model 
fits (solid curves) (see text). 
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4. Non-Identical Two-Particle Short-Range Correlations 

Lednicky et al. jL0) pointed out that measurements of short-range relative momen- 
tum correlations in the co-moving pair frame of non-identical particles could be used to 
constrain dynamical models of the emission source, in particular the average space-time 
differences between emission regions for different particle species. f\ For example, for ir + K~, 
the attractive Coulomb interaction produces correlations greater than unity. Correlations 
obtained for pairs where the pion speed is selected to be greater than the kaon's are de- 
noted and vice versa for C_. If the pion and kaon emission times are the same, 
for example, observation of C+ K /CZ K > 1 implies that the average ir + emission radius is 
smaller than that of the K~ and vice versa if C1 K /C^_ K < 1. Similar arguments apply to 
the other charge sign combinations. 

Three-dimensional Pratt-Bertsch |1 tt K correlation ratios for the top 12% central 
events from the 130 GeV Au+Au collision data from STAR were measured and analyzed 
(see contribution to these proceedings by F. Retiere). Transverse velocity magnitudes 
(/?) were restricted to the range from ~ 0.7 — 0.75c due to limitations imposed by dE/dx 
based particle identification with the TPC. Due to the inherent ambiguity in the analysis 
the data are consistent with pion emission which in one limit occurs at relatively smaller 
transverse position (~ 4.1 fm in the lab frame) than the kaons, assuming equal emission 
times, or in another limit occurs an average of 5.4 fm/c later than kaons assuming emission 
at equal transverse positions. 

These 7r K correlation ratios can be understood in terms of the "blast-wave" model 
discussed above. The selection of /3 ~ 0.72c particles together with the brief emission 
duration and sharp radial cut-off of the emission source geometry in the blast-wave model 
results in the emission zone being a thin shell near the outer edge of the source. Thermal 
smearing effects, which are larger for lighter mass particles (pions), conspire with the 
latter effect to cause the average pion emission radius to be smaller than the kaon's. The 
parameters of the blast-wave model which describe these data [T = 110 MeV, (p) = 0.6 
(average transverse flow rapidity), R = 13 fm, At = 1.5 fm/c and t$ ph 10 fm/c], are 
consistent with those which describe the pt distribution shapes, pion HBT radii, and 
R /R s ratio. 

5. Long-Range Correlations in Pseudorapidity Space 

Two-particle momentum space correlations in general are six-dimensional objects. De- 
pending on the nature of the dynamics to be studied and possible limitations of statis- 
tics, projections onto fewer dimensions, usually one or two, are typical in most analy- 
ses. For the STAR data presented, charge independent (CI) and charge dependent (CD) 
combinations of charged-sign pairs were formed. The CI correlations were obtained by 
summing over all pairs regardless of charge sign. The CD correlations are defined by 
C(CD) = {C ++ + C ) — (C+_ + C_ + ), where C±± are charge-pair specific correlations. 

Correlations in pseudorapidity (77) space are affected by string fragmentation, source ex- 
pansion, medium effects, etc. Fig. |3| displays, for example, the two-particle CD pseudora- 
pidity correlation rji<S>T]2, for 130 GeV Au+Au collision data from STAR for central events 



3 This analysis method cannot simultaneously determine spatial and temporal differences. 
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(multiplicities > 0.81 of maximum) using primary particles with 0.15 < px < 2 GeV/c, 
| T] | < 1.3, and full azimuthal acceptance. The charge-pair specific correlations were ob- 
tained by dividing the number of particle pairs of each charge combination from the same 
event (sibling pairs) in a given bin (0.1 in width), summed over all events in the sample, 
by the number of corresponding charged particle mixed-pairs in the bin from different 
events. The latter events were required to have primary vertex locations within 5 cm in 
the TPC along the beam axis. The correlation function was normalized based on the total 
number of pairs used in the numerator and denominator. Large-scale structure covering 
the full STAR pseudorapidity acceptance is evident in Fig. |3| where the dark band along 
the rji = r]2 diagonal is negative and the off-diagonal corners are positive. The dominant 
structure is with respect to the difference variable 771—772- Projection onto the 771—772 axis, 
shown in Fig. f|, results in the CD autocorrelation. This is one example of long-range axial 
(longitudinal or azimuthal) correlations which potentially carry new information about 
the dynamics and evolution of the collision. 
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Figure 3. Two-dimensional r\\ ® r] 2 charge 
dependent correlations for ^/snn = 130 GeV 
Au+Au collisions. 
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Figure 4. CD autocorrelation projec- 
tion onto the (r)i — r] 2 ) axis for ^/snn = 
130 GeV Au+Au collisions from Fig. [| 



6. Long-Range Azimuthal Angular Correlations 

Two-dimensional azimuthal angle correlations, 0i <8> 02, are dominated by elliptic flow 
contributions which produce a l + 2u| cos[2(0i — 2 )] dependence. Autocorrelation projec- 
tions onto the <f>% — 02 axis are shown in Figs. |5] and ^| for like and unlike charge-sign pairs 
for ^snn = 130 GeV Au+Au collision data from STAR for combined mid-central and 
central events (multiplicities greater than 0.47 of maximum) assuming the same primary 
track acceptances described in the preceding section. In addition to elliptic flow, this kind 
of analysis reveals long-range non-flow components [|ll| (e.g. non-collective, few-body pro- 
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Figure 5. Autocorrelation projection onto Figure 6. Autocorrelation projection onto 

the (0i — 02 ) axis i n radians for like-sign the (0i — 2 ) axis in radians for unlike- 

charged particle pairs for y/s NN = 130 GeV sign charged particle pairs for y/s NN = 

Au+Au collisions. 130 GeV Au+Au collisions. 



cesses such as resonance decays, mini-jets, etc.) which appear as the nan — cos[2(0x — 02)] 
dependences in Figs. [5] and |J 

A possible method to assess the importance of non-flow components has been discussed 
by Borghini et al. ||12|| . It is based on the technique of forming four-particle cumulants 
for studying azimuthal anisotropy in particle production from heavy ion collisions. In the 
cumulant expansion, dependences on lower-order correlation contributions are eliminated 
term-by-term leaving only those contributions arising from true, four-particle correlations. 
In this way, the deduced values of v 2 obtained with this procedure are expected to more 
accurately measure the collective, elliptic flow produced in heavy ion collisions. 

This new method was implemented in STAR JTIJ and used to analyze the ^s^n = 130 
and 200 GeV Au+Au data. Values of v 2 obtained in this way are reduced about 15% 
relative to that obtained with either the standard reaction plane method ]1^] or with two- 



particle angular correlations [II|. A comparison of v 2 versus centrality using four-particle 
cumulants for the ^/snn = 130 and 200 GeV Au+Au data is shown in Fig. ^ where it is 
noted that v 2 on average increases about 7-10% at the higher energy. 

From the study of elliptic flow in STAR it has been found that the event plane resolu- 
tion is sufficient to conduct three-dimensional Pratt-Bertsch HBT analysis with respect 
to event plane. This type of analysis of charged pion data for ^/s^n = 130 and 200 GeV 
Au+Au collisions was reported in these proceedings by M. Lopez-Noriega. This analysis 
indicates an out-of-event-plane extended source with ~ 5% azimuthal anisotropy for the 
geometrical source shape. The azimuthal dependence of the HBT radii cannot be ex- 
plained assuming only azimuthal momentum anisotropy or a source with only azimuthal 
anisotropy in the density In particular, the blast-wave model provides consistent de- 
scriptions of both v 2 (m,pT) and the azimuthally dependent HBT radii. Interestingly, the 



azimuthal source shape anisotropics at RHIC and the AGS [14] are both out-of-plane 



extended but the magnitude at RHIC is about one-tenth as large. 
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Figure 7. Preliminary azimuthal anisotropy 
measurements of t>2 versus centrality for 
nonidentified charged particles using four- 
particle cumulants for ^/snn = 130 and 
200 GeV Au+Au collisions. 



Figure 8. Perspective view of the two- 
dimensional X(mTi) <8> X(rriT2) charge 
independent correlations for ^snn = 
130 GeV Au+Au collision data. 



7. Long-Range Transverse Momentum Correlations 

The shape of the transverse momentum spectrum of charged hadrons near mid-rapidity 
produced in relativistic heavy ion collisions is affected by quantum statistics, thermaliza- 
tion, transverse expansion, hard QCD processes, mini-jets, etc. In addition, non-statistical 
or dynamical fluctuations within each event or from event-to-event modify the shape of 
the inclusive pt distribution, increasing the yield at higher px- For example, an intrin- 
sic Boltzmann distribution with gaussian-like event-by-event fluctuations in the inverse 



temperature parameter becomes a Levy distribution in the inclusive pr spectrum JL5 
Similarly, fluctuations in inverse temperature produce a long-range "saddle shape" in the 
two-dimensional m^i ® correlation function [15| , which is the transverse momentum 



space analog of the two-dimensional correlations shown in the preceding sections. This 
is shown in Fig. |] for the 5% most central Au+Au events at a/saT/v = 130 GeV using all 
charged particles with 0.15 < px < 2.0 GeV/c, \t)\ < 1.3, and full 2tt acceptance.^ The 
dominant features in these data include HBT and Coulomb interaction correlations which 
produce the ridge at low %, the peak at large momentum, plus the long-range saddle 
shape. All features have similar amplitudes and are the result of non-statistical processes. 

Further evidence that non-statistical event-by-event fluctuations occur in RHIC colli- 
sions is observed in the distribution of event- wise average pt ((Pt))- From application of 
the central limit theorem (CLT) |i5|Jl6[| it is known that in the absence of non-statistical 



4 Approximately constant statistics in each bin were achieved by mapping px into X(m<r) where X(my) = 
1 — [(1 + mx/T) exp(— tut /T)]/[{1 + wlq/T) exp(— mo/T)]; ttlq = pion mass was assumed. 
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fluctuations, the mean and rms width of the event-wise (py) distribution are equal to pPf 
and cjp T /\/ r N, respectively, where pr and <t| t are the mean and variance of the inclusive 
Pt distribution and N is the event multiplicity. Increased width in the data beyond this 
amount provides further evidence for non-statistical fluctuations. 

The inclusive pt spectrum of charged particles near mid-rapidity over the pr range of 
interest here is described very well by a Gamma distribution JI6J. The (pt) distribution 
for RHIC events which would be expected in the absence of non-statistical fluctuations is 
accurately described by the iV-fold convolution of the inclusive px Gamma distribution. 



This is also a Gamma distribution [16j whose variance is in agreement with the CLT |T5 
16] . The accuracy of the iV-fold Gamma distribution for representing the statistical- 
fluctuations-only (pt) reference distribution is assured, for moderate to large values of N, 
by the strong suppression of higher-order moments in this distribution. This suppression 
is driven by powers of N which increase with moment number. Traditional mixed event 
and Monte Carlo sampled references were also obtained and shown to be consistent with 
the more accurate Gamma distribution reference. 

The frequency distribution^ of event-wise (pr) for the upper 15% central Au+Au colli- 
sions at y/SNN = 130 GeV from STAR is shown in Fig. |9] in comparison with the Gamma 
distribution reference. These data are from 183K events and use 70% of all primary 
charged particles (less than 100% due to tracking inefficiency and track quality require- 
ments) with 0.1 < pt < 2.0 GeV/c, \r]\ < 1, and full 2n acceptance. Striking evidence is 
observed for large, non-statistical fluctuations which increase the width of the measured 
distribution by 14±0.5(stat)% relative to the Gamma distribution reference. 

Quantitative measures of (pt) non-statistical fluctuations are provided by a number of 
quantities [|TJ,|T^,^], for example: 

1 e 

Aa l T = tE^j «Pr)i - Pt? - 4 T = 2a PT Aa PT , (3) 

1/2 

~°PT, (4) 



1 1 



°"<p T >, dynamical = ~Yj _ T) ^ (S I>I, *I>I ,- (5) 

j=l j\ j J i^Li'=\ 

where e is the number of events, j is the event index, Nj is the event multiplicity, (N) is 
the mean multiplicity, i is a particle index, and 5px , = Pt t — Pt- For minimal variations 



of Nj within the event ensemble 



Aa pT — $ PT — — a (pr), dynamical- (6) 

Charge sum (CI) and difference (CD) measures corresponding to Aa PT are given by 

(N)^ 2 Pt ,ci/cd = (N + )Aa 2 p T++ + (M> Aa 2 pT __ ± 2y/(N + ) (N_)Aa 2 pT+ _, (7) 
5 The distribution shown is for the quantity [(pt) — PtI/Wpt/^N] which is designed to suppress the 



slight dependence of px on N, remove the l/y/N statistical broadening, and eliminate dependence on 
the variance of the inclusive distribution. 
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Figure 10. Mean p? CI (solid triangles) and 
CD (open triangles) fluctuation measures for 
130 GeV minimum bias Au+Au events versus 
centrality as explained in the text. The CD 
quantities have been multiplied by 3 for clar- 
ity. 



where AcTp r±± = 2ap T ±Aa PT ±± 1 A(Tp T+ _ = 2^/a PT+ a PT ^Aa PT+ ^, and AcTp T+ _ is evaluated 
noting that the covariance, cr| T+ ^ T _, vanishes via the CLT. 

The centrality dependence of Aa PT is shown in Fig. [H] for 205K ^/s^n = 130 GeV 
Au+Au minimum bias events from STAR using 70% of all charged primary particles for 
charge independent and charge dependent measures. Statistical errors are ±0.5 MeV/c. 
The dashed lines are polynomial fits to the data, the solid lines are estimated extrapola- 
tions to 100% of the primary particles in the acceptance, and the shaded bands indicate 
±15% systematic errors. 

These data provide a quantitative measure of the non-statistical fluctuations visually 
demonstrated in Fig. |^ and reveal an intriguing dependence on centrality. They also 
show that the CI (CD) fluctuations at RHIC are larger (smaller) than at the SPS p0 



The PHENIX experiment's null result [21], which was carried out with much smaller 
rj, (p acceptance than in STAR, does not contradict STAR's observation of strong, non- 
statistical fluctuations. 



8. Summary and Conclusions 

During the first two years of RHIC operation the STAR experiment has produced 
a wealth of new correlation, flow, and event-by-event fluctuation measurements.^ The 
results for the ^/snn = 200 GeV Au+Au collisions are similar to those for the 130 GeV 

6 Event-by-event fluctuations in net charge were also measured for the 130 and 200 GeV Au+Au and p+p 
collision data; see the contribution to these proceedings by C. Pruneau. 
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data. Compared to heavy ion results at lower energies from the AGS and SPS, the 
RHIC events display (1) similar HBT correlation lengths, (2) R Q /R S « 1, (3) larger pion 
phase space densities, (4) stronger elliptic flow t> 2 , (5) larger (smaller) charge independent 
(dependent) mean fluctuations, (6) strong transverse expansion with (apparent) brief 
emission duration, and (7) similar out-of-event-plane extended sources but with much 
smaller azimuthal source shape anisotropies. The STAR measurements, the successful 
description of the data in terms of the blast-wave model, and the general success of 



the hydrodynamical approaches |22[ at RHIC energies, taken together, suggest that the 
Au+Au collisions at RHIC are governed by rapid thermal equilibration, very high initial 
pressures, followed by hydrodynamical expansion. Analyses of the emerging fluctuation 
measurements from STAR will also help determine the source structure and dynamics 
and allow further differentiation among QGP and purely hadronic scenarios. 
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